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Abstract. We use a semiclassical three-fluid model to analyze the conditions for 
spatial phase separation in a mixture of fermionic 6 Li and a (stable) Bose-Einstein 
condensate of 7 Li atoms under cylindrical harmonic confinement, both at zero and 
finite temperature. We show that with the parameters of the Paris experiment [F. 
Schreck et ai, Phys. Rev. Lett. 87 080403 (2001)] an increase of the boson-fermion 
scattering length by a factor five would be sufficient to enter the phase-separated 
regime. We give examples of configurations for the density profiles in phase separation 
and estimate that the transition should persist at temperatures typical of current 
experiments. For higher values of the boson-fermion coupling we also find a new phase 
separation between the fermions and the bosonic thermal cloud at finite temperature. 
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1. Introduction 



Following the achievement of Bose-Einstein condensation in alkali gases |]T|, a great 
deal of experimental effort is being devoted to trap and cool gases of fermionic isotopes 
of alkali atoms. The usual evaporative cooling scheme is ineffective, since the Pauli 
exclusion principle forbids s-wave collisions between spin-polarized fermions. This 
problem has been circumvented by recourse to sympathetic cooling, which relies on 
s-wave collisions with atoms belonging to different species, and several groups have 
thus reached in the Fermi cloud temperatures as low as T = 0.2 T F P, [| |J, T F being 
the Fermi temperature. In particular, an experiment in Paris has produced a mixture 
of a degenerate Fermi gas and a stable Bose-Einstein condensate, made from Lithium 
isotopes ||. 

The condensate-fermion mixture is an interesting system for the study of the effects 
of the interactions, since it can show spatial phase separation of the two components 
on increasing the value of the boson-boson and boson-fermion coupling constants ||. 
This is an example of a quantum phase transition ||. Several configurations for phase- 
separated clouds have been proposed in . A schematic phase diagram for the trapped 
boson-fermion mixture at zero temperature has been given in ||, while the phases and 
the stability of a homogeneous boson-fermion mixture have been studied in 0. In the 
experiments it might be possible to investigate the phase-separated region by exploiting 
optically or magnetically induced Feschbach resonances to increase the values of the 
boson-boson and/or the fermion-boson scattering length fTT)|j . 

In this letter we have chosen the trap geometry and parameters typical of the 
Paris experiment and have analyzed the possible phase-separated configurations; we 
also indicate how various equilibrium observables are affected by phase separation. 

The paper is organized as follows: in Sec. ||we present the model which we have used 
and discuss its limits of applicability; Sec. |] gives the conditions for phase separation at 
zero temperature and illustrates the column- density profiles corresponding to different 
phase-separated configurations. The results are extended to finite temperature in Sec. |. 



2. The model 



We adopt a three-fluid model for a mixture of bosons and fermions at finite temperature 
T [PH . The model is used to evaluate the density profile of the condensate (n c ) within 
the Thomas- Fermi approximation and those of the fermionic (n/) and thermal bosonic 
component (n nc ) within the Hartree-Fock approximation. 

We shall assume that the number of bosons in the trap is large enough that the 
kinetic energy term in the Gross-Pitaevskii equation for the wave function ?/>(r) of the 
condensate can be neglected fll2| . The density profile of the condensed atoms is 

n c (r) = ^ 2 (r) = - V h ex \v) - 2gn nc (r) ~ fn f {v)]/g (1) 

for positive values of the function in the brackets and zero otherwise. Here, g = 
Anffabb/mb and / = 2nh 2 abf/m r with a&& and a^f the boson-boson and boson-fermion s- 
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wave scattering lengths and m r = m^mf / (m^ + mf) with m& and nif the atomic masses; 
\ib is the chemical potential of the bosons. The system is confined by external axially 
symmetric potentials 

V h f{v)=m bJ ul f {r> + \l f z*)/2 (2) 

where u>bj are the frequencies and Xbj the anisotropics of the traps. 

As already proposed in early work on the confined Bose fluid [0, we treat both the 

/efj 
b,f 



thermal bosons and the fermions as ideal gases in effective potentials V bf (r) involving 



the relevant interactions. We write 

V b eff (r) = V b ext (r) + 2gn c {v) + 2gn nc {v) + fn f (r) (3) 

and 

Vf{v) = Vf{v) + fn c (r) + fn nc (r). (4) 

We are taking the fermionic component as a dilute, spin-polarized Fermi gas: the 
fermion-fermion interactions are then associated at leading order with p-wave scattering 
and are negligible at the temperatures of present interest |14| . Here and in Eq. ([!]) , the 
factors 2 arise from exchange. 

We may then evaluate the thermal averages by means of standard Bose-Einstein and 
Fermi-Dirac distributions, taking the thermal bosons to be in thermal equilibrium with 
the condensate at the same chemical potential and the fermions at chemical potential 
fif. In this semiclassical approximation the densities of the cold- atom clouds are 

d 3 P ' " •' ■ 1 



The chemical potentials are determined from the total numbers of bosons and fermions. 

The above model is valid for bosons when the diluteness condition n c a\ b 1 
holds. For the Fermi cloud in the mixed regime the condition fc/Ob/ 1 with 
kj = (A&Np) 1 ^ / aho should hold, but this is not a constraint in the regime of phase 
separation where the boson-fermion interaction energy drops rapidly (see below). 



3. Conditions for phase separation 

Phase separation is complete when the overlap between the bosonic and the fermionic 
density profiles vanishes. The transition to this regime depends on the confinement and 
on the numbers, masses and interaction strengths of the trapped fermions and bosons. 
In the case of equal numbers of bosons and fermions and equal trapping frequencies the 
condition of phase separation can be written as 

37r fa \ ^ 

kfabb> iw f ) (6) 

at zero temperature |8], . 

Let us focus on the geometry and on typical values of the numbers of atoms in the 
Paris experiment || for the 7 Li- 6 Li mixture: u b = 2n x 4000 Hz, ojf = 2n x 3520 Hz, 
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Figure 1. Boson-fermion interaction energy (in units of Nfruib with N = 20000 being 
the total number of particles), as a function of \n(a,bf / abb)- The curves are obtained 
for abb = 5.1 ao and equal numbers of bosons and fermions at various values of the 
temperature: T = (solid line), T = 0.2 Tp (dashed line) and T = 0.6 Tp (dotted- 
dashed line, on a scale reduced by a factor 20). The arrow indicates the location of 
the Paris experiment. 



X b ~ Xf ~ 1/60, and N b ~ Nf ~ 10000. Since the two trap frequencies are of 
the same order of magnitude, we can use Eq. (|(J) to estimate at which values of the 
scattering lengths the phase transition occurs. For the above experimental parameters 
the condition (g) becomes £ = alj/a bb > 7000 a Q . Using the predicted values for the 
scattering lengths of the 7 Li- 6 Li mixture, that is a bb = 5.1 ao and aj,/ = 38 ao @, we 
have £ ~ 283 a and the gases are not as yet in the phase-separated regime. However, 
in a Feschbach-resonance experiment the scattering lengths could be increased. There 
is therefore the possibility to fulfill Eq. @ and to enter the phase-separated regime. 

Since the positions of the Feschbach resonances for a b f and a bb would not necessarily 
coincide, in order to explore the phase diagram for the boson-fermion mixture we have 
first of all kept a bb fixed at its non-resonant value and changed a b f. In this case, in 
order to reach phase separation it is sufficient to increase the value of the boson-fermion 
scattering length by a factor five, namely abf > a b f ~ 200 ao- 

To check this estimate we have evaluated the boson-fermion interaction energy 

Eint = f J d 3 r [n c (r) + n nc (r)]n/(r) , (7) 

using the model described in Sec. 0. The behaviour of this quantity at T = as a 
function of a b f/a bb , with a bb = 5.1ao, is shown in Fig. [l| (solid line). Because of finite 
size effects the transition to the phase-separated regime shows a large crossover region 
and the interaction energy goes smoothly to zero when separation is complete [15| . 
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Figure 2. The symmetric configuration at the threshold of the phase-separation 
regime, obtained for = 5.1 do and a>bf = 200 a . Top panels: column densities in 
the {x, z) plane of the fermionic cloud (left) and of the condensate (right); bottom 
panels: topview images of the same profiles. The size of the figures in the {x, z} plane 
is 10/im x 665/im. 




Figure 3. The fully asymmetric configuration for abb = 5.1 ao and a&/ = 1000 ao- Top 
panels: column densities in the {x, z] plane of the fermionic cloud (on the left) and 
of the condensate (on the right); bottom panels: topview images of the same profiles. 
The size of the figures in the {x, z} plane is 10/im x 665/im. 
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For the value a&/ = 38ao corresponding to the current experiments (marked by the 
arrow in Fig. [[]), the boson-fermion interaction energy is close to its maximum value. 
We find that the overlap between the two clouds vanishes if a^j > 200ao, in accord with 
the estimate given just above. 

To make close contact with experiments we have calculated the column densities of 
the fermionic and condensate vapours in the phase-separated regime for various values 
of a^f. At the threshold value we have found only one possible configuration, with 
a central hole in the fermionic cloud which is occupied by the condensed bosons (see 
Fig. For this configuration there is only a quantitative difference in the value of the 
contrast at the center of the fermionic cloud as compared to the mixed state. 

The regime of phase separation is univocally characterized by the presence of other 
configurations at higher values of a&/. Starting from different initial conditions, for 
1000 ao < abf < 2000 a we have found, in addition to a configuration which is similar 
to that in Fig. 0, another with a different ellipticity of the hole and a third one which 
is asymmetric in the axial direction (see Fig. ^). For the values of the parameters that 
we have used the lowest energy configuration is the first (symmetric) one. However, 
also the other metastable states are close in energy and could perhaps be realized in 
the laboratory by applying an adiabatic perturbation having the symmetry of the final 
state. The asymmetric state could easily be identified in the column-density pictures. 

Other interesting configurations would be observed in an experiment in which both 
scattering lengths could be changed. On increasing a bb the size of the condensate 
increases towards that of the fermion cloud and the hole in the latter widens, until 
a situation is reached in which the fermions are pushed out in the radial direction and 
make a ring around a cigar of bosons. Another peculiar shape is a "double sandwich" of 
bosons, in which two slices of condensate are separated by a central slice of fermions and 
surrounded by a fermionic shell. As a last remark, the energy of all these configurations 
depends on the values of the scattering lengths: it may thus happen that the most stable 
state is the asymmetric one or has other exotic shapes [ITH . 



4. Effects of temperature 

At low temperature (T ~ 0.2 Tp) the density profiles for the various configurations are 
essentially the same as those found at T = 0, with the addition of a small bosonic 
thermal cloud. For values of the boson-fermion scattering length around the critical 
value abf the phase separation involves only the condensate and the fermionic cloud: 
the total boson-fermion interaction energy E int has a minimum but does not vanish 
(see the dashed line in Fig. [I]), since the fermionic and the bosonic thermal clouds still 
overlap. However, with a further increase in the scattering length, up to values of order 
10 5 ao, the model yields full phase separation between fermions and bosons. 

At higher temperature (T ~ 0.6 Tp) we do not find phase separation between the 
fermionic cloud and the condensate (see the dotted-dashed line in Fig. p. Separation 
between the total bosonic gas and the fermions is instead found in this model at 
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Figure 4. Variance of the axial distribution of a cloud of 10000 atoms of 6 Li taken 
at x = (in units of Zj = 2k B T F /m f uj 2 f ) against T/T F with T F = {<oN f ) l ^fva) f /k B . 
The experimental data (filled circles; from Ref Q and private communication) are 
compared with the behaviour of the classical gas (dotted line) , of the ideal Fermi gas 
(full line), and of the fermionic cloud interacting with its bosonic partner (Nb = 20000) 
using the values = 5.1 ao and abf — 200 do (dotted-dashed line). 



a,bf ~ 10 5 ao- In this case a symmetric configuration is generated, in which the fermions 
are pushed outside the whole boson gas. The size of the hole in the fermionic cloud is 
determined by the width of the bosonic thermal cloud, making the finite-temperature 
transition much more evident in the column-density images than at T = |I6| . 

Another quantity which can be measured as a function of temperature is the width 
of the fermionic density profile in a given direction, e.g. in the z direction (see Fig. ^). 
While for the values of the scattering lengths in the Lithium experiment the width is 
predicted to be almost the same as that of the non-interacting Fermi gas [XT], larger 
increases in width are found in some phase-separated configurations. This observable 
is particularly sensitive to the configurations which involve a large displacement of the 
fermions, such as the "double boson sandwich" 
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In conclusion, we have examined the main observable features of spatial phase 
separation in gaseous mixtures of Li isotopes, with specific focus on the setup of the 
Paris experiments. A full account of our results will be given elsewhere [[L6| . 

We acknowledge support from MIUR under PRIN2000 and by INFM under 
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